Bead on a rotating circular hoop: a simple yet feature-rich 

dynamical system 

Shovan Dutta^ and Subhankar Ray^ 
^ Department of Electronics and Telecommunication Engineering, 
Jadavpur University, Calcutta 700 032, India and 
^ Department of Physics, Jadavpur University, Calcutta 700 032, India 

The motion of a bead on a rotating circular hoop is investigated using elementary 
calculus and simple symmetry arguments. The peculiar trajectories of the bead at 
different speeds of rotation of the hoop are presented. Phase portraits and nature of 
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study introduces several important aspects of nonlinear dynamics. It is suitable for 
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I. INTRODUCTION 



A bead moving on a rotating circular hoop is a classic example studied in several text- 
books on classical meclianic^^ and nonlinear dynamics. It exhibits various modes of mo- 
tion, including some peculiar ones, such as, oscillations confined to one side of the hoop 
and complete revolutions, for appropriate initial conditions. It also shows a wide array of 
features of dynamical systems. It is useful for demonstrating different classes of fixed points, 
bifurcations, reversibility, symmetry breaking, critical slowing down, trapping regions, ho- 
moclinic and heteroclinic orbits and Lyapunov functions. It is also an interesting example of 
a constrained system, illustrating the use of Lagrange multipliers for determining constraint 
forces. 



In this article we investigate in detail, the dynamics of this bead-hoop system. The 
number and nature of equilibrium points alter with change in speed of rotation of the hoop. 
This results in some extraordinary modes of motion, phase trajectories and bifurcations. At 
a critical speed of rotation of the hoop, an interesting relation between the time period and 
amplitude of oscillation of the bead is observed. 



Our study is based on elementary calculus and simple symmetry arguments. A student 
reader with a background of classical mechanics and basic calculus will find the study both 
accessible and interesting. 



The article is organized as follows. In section |TT| the physical system is described. The 



equation of motion of the system is derived from the Lagrangian. In section III , the effective 
potential is analyzed to find the equilibrium positions of the bead for all possible speeds of 
rotation of the hoop. The distinct modes of motion of the bead for different initial conditions 
are described and numerical plots of its trajectories are presented. The constraint forces are 
determined from the method of Lagrange multipliers in section |IV[ In section |V} the nature 
of the fixed points of the system are analyzed from symmetry properties of the system. The 
phase trajectories and bifurcations are examined. Finally, we conclude with a discussion on 



connections to other systems in section VI 
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II. THE PHYSICAL SYSTEM 

A bead of mass m, moves without friction on a circular hoop of radius a. The hoop 
rotates about its vertical diameter with a constant angular velocity uj. The position of the 
bead on the hoop is given by angle 9, measured from the vertically downward direction {—z 
axis), and cj) is the angular displacement of the hoop from its initial position on the a;-axis ( 
Figure [T]) . 




FIG. 1. Schematic diagram for bead on the hoop 



The kinetic and potential energies of the bead are given by, 



T 



1 



-ma 



2(n2 



+ sin 



2 /D J,2\ 



and V = —mga cos 9 



(1) 



respectively, where g is the magnitude of the acceleration due to gravity. The Lagrangian 
of the system is. 



ma 



+ sin^ 9) + mga cos 9 , where uj 



(a constant). 



Using the Euler-Lagrange equation, the equation of motion is obtained as, 

9 = sin6'(c<;^ cos6' — g/a). 



(2) 



(3) 



Denoting = g/a, k = uj'^/uj'^ and defining t = Uct, the equation of motion may be written 
in dimensionless form as, 

d^9 



9 = sm9{kcos9 — 1) , where 6' 



(4) 



The symmetry of the system about the vertical axis manifests in the invariance of the equa- 
tion of motion when 9 is changed to —9. The system also exhibits time reversal symmetry 
as ^ is invariant under the transformation r — ?■ — r. 
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Let us denote v = 6 and write \4m as 



which upon integration yields, 



dv 

V — = — sin 6* + sin0 cos^, (5) 
dO 



1 k 

-t,2 = cos ^ - - cos 29 + A, (6) 

where A is a constant of integration. We can identify a conserved quantity, e as, 

e= ^fcsin^^-cos^ . (7) 

It may be termed as the effective energy that corresponds to the system described by (|4]). 
However, it is to be noted that, the mechanical energy of the original bead- hoop system 
is not conserved due to the work done in preserving the constant angular velocity of the 
hoop. Equation ([T]) provides a natural separation of the effective energy into effective kinetic 
energy 16*'^ and effective potential energy, U{6), given by, 

u{e) = -hisixi^e -cose . (8) 

As U{6) is an even function of 6, we restrict our attention to G [0,7r]. O'"^ is completely 
determined by 6 and e. The sign of 6\ however, cannot be determined uniquely from ([t]). 
The parameter k grows with increasing angular velocity u of the hoop and A; = implies 
that the hoop is stationary. 

It is useful to consider the first and second derivatives of U{9), 

U'(9) = ^ = -sinOik cos 9-1) , and (9) 
do 

U"(9) = ^ = -(kcos2e -cose) . (10) 
do^ 

Comparison of ^ and ^ shows that, 

e" = -u'{e) = -sme{i-kcose). (ii) 

The effective potential has two extrema at 6 = and 6 = n, where the first derivative of the 
potential vanishes. Evaluating the second derivative at these points, one finds that = is 
a minimum {U (0) > 0) and e = TT is a maximum {U (0) < 0) of the effective potential as 
shown in Figure |2} Due to the right-left symmetry of the hoop, one concludes, 6 = —tt is 
also a maximum point of the effective potential. At these points the values of the effective 
potential are U{0) = Umin = — 1 and f/(±7r) = Umax = 1- 
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FIG. 2. Minimum of U{9) is flatter and maxima sharper as fc — >• 1. The soUd curve is for k = 0.010, 
dashed curve is for k = 0.998 

III. MODES OF MOTION AND TRAJECTORIES OF THE BEAD 

A. < k < 1 

When k = 0, the effective potential takes the form of the potential energy of a pendulum. 
As k increases from zero towards 1, U" (0) — )■ and |t/"(±7r)| — )■ 2, implying that the 
minimum at 6' = is flatter and the maxima at ^ = ±7r are sharper (Figure |2]). 

For k ^ 0, at different effective energy values, the actual motion of the bead will depend 
on the initial conditions which we discuss below. 
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FIG. 3. Effective potential U{e) versus 6* for < /c < 1 
(i) Effective energy, e = —1 
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^'^ = 2(e — U{6)) must be positive, hence the lower bound of e is equal to the minimum 
value of U{6) which is —1 (Figure Isl). This corresponds to the initial condition ^(0) =0 and 



^'(0) = 0. From (11) however, 9" = as well, independent of the value of k. Hence, the 
bead stays at the bottom of the hoop {9 = 0) for all times. 



(ii) -1 < £ < 1 

The angular position 9 of the bead oscillates between —6*0 and +6*0, corresponding to the 
turning points B and C shown in Figure [sj At these points 9' = 0, and the effective potential 
equals the effective energy. From (14|, 0" < and this restoring acceleration ensures that 9 



ff 




FIG. 4. Out of phase oscillation of O' and O" for k = 0.75, = 27r/3 

approaches 0. hi 9 = 0, acceleration 9" is zero, whereas the velocity 9' attains its maximum 
magnitude. Hence, the bead slides past the bottom of the hoop and 9 decreases further, 
becoming negative. When < 0, 9" becomes positive, and this acceleration reduces the 
magnitude of 9\ making it zero at —9q. 9" remains nonzero and positive at —9^. The 
acceleration 9" is positive for all negative values of 9. Kt 9 = —9q, 9" > and velocity is 
momentarily zero. From the next instant, velocity begins to increase as 9" > 0. Thus 9' and 
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9 oscillate out of phase with each other as shown in Figure ^ At -^0, ^ > 0, ^ becomes 
positive, 9 begins to increase. The bead retraces its path reaching 9q and the process is 
repeated. The amplitude of oscillation increases with e. The time period of the oscillation 
of the bead between ±6*0 is given by, 

"^'^ d9 



T 



(12) 



'0 a/(cos6' — cos6'o)(2 — k{cos9 + cos 6*0)) 
Figure [5] shows the trajectory of the bead after the elapse of different time intervals, with 
/c = 3/4 and amplitude 9q = 2tt/3. As the bead goes through half an oscillation (with 
amplitude 27r/3 and k = 3/4), the hoop rotates through an angle given by, 

[a T\ \/k^ 
^2 = 



A(j) = co 



-T = 3.539 radians 




(a) r = 28 (b) r = 62.5 (c) r = 62.5, top view 

FIG. 5. Trajectory for k=0.75 , 6(0) = 27r/3, 6' (0) = at different r values 



(iii) e = 1 

At = ±7r, the effective potential energy is also +1. Thus the effective kinetic energy and 
9' is zero. From the expression of acceleration in (|4]), we find = at = ±tt. This agrees 
with the fact that 9 = ±tt are the position of maximum potential energy. From ([T]), we 
obtain, 

^" =4cos2^[l + A;sin2^] . (13) 

The time taken to reach the highest point [9 = vr) from an initial point ^0 may be computed 
as, 



T 



1 r 



d9 



^ Jeo cos 



(14) 
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The integral can be transformed with the substitution u = n — 6, as, 

1 r"^° du , , 
T = - = 15 

2 7o sin (m/2)v/1 + A;cos2 (m/2) 

In the neighbourhood of m = 0, corresponding to 6 close to vr, the integrand behaves as 
and hence, the integral diverges. This implies that the bead approaches 6 = tt 



infinitesimally slowly and takes an infinite amount of time to reach the top of the hoop 
(Figure IgI). By virtue of time-reversal symmetry and the left-right symmetry of the hoop, 



FIG. 6. Trajectory for k = 0.75 , 9(0) = 27r/3 and e = 1 

for ^'(0) < 0, the bead reaches the bottom where 9' has the maximum magnitude of 2 (from 
([7j)) and then asymptotically approaches 6 = — vr at ever-decreasing speed. 

(iv) e > 1 

As the maximum value of U{6) = 1, e — U{6) > and hence 6' is not zero for any 6. Thus 
if the bead starts initially with ^'(0) > and at a position ^(0) G (0, vr), then ^ "(0) < and 
9' starts to decrease. 9' attains its minimum value when the bead reaches 9 = n. At this 
position, 9" = 0, however, 9' is not zero and 9 continues to increase. Hence the bead moves 
past the top of the loop. For 9 G (vr, ^), 9" is positive and 9' increases, becoming maximum 
when the bead reaches the bottom of the hoop. Thus 9 increases or decreases indefinitely, 
depending on the sign of 9' (0) and the bead slides over the entire hoop periodically (Figure 
Q. For ^'(0) < 0, the same rotation occurs in the opposite sense. The maximum and 
minimum values of 9' can be calculated from the expression of the effective energy. The 
time period of rotation is given by. 
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When |6''(0)| is very large, nearly all the energy of the bead is the kinetic energy of rotation. 
Hence the influence of gravity becomes small in comparison, and this rotation may be 
considered to be uniform. In this limit, the time period approaches the value 27r/|^^'(0)|. 




(a) r = 20 (b) r = 40 (c) r = 200 

FIG. 7. Trajectory for k=0.01, 6^(0) = 27r/3, 9' {0) = 1.4 for different r. 



B. k = 1 

U{9) is maximum at ^ = vr as before. However, at ^ = 0, both U' (0) and U" (0) vanish. 
On expanding U{9) about 9 = 0, we get, 

U{9) = -1 + ^ + 0(^^) . 

For 1^1 > 0, in the neighbourhood of 9 = 0, U{9) is always greater than f/(0). Hence ^ = is 
a minimum of the potential. The potential is more fiat than for the case k < 1. Depending 
on the values of e, the bead can undergo all the types of motion described in the previous 
subsection. However, for oscillatory motion, the period of oscillation for a given amplitude 
is larger when k = 1 than when k < 1. This becomes clear if ^ is expanded in the vicinity 
of ^ = 0, 

^" = (^-f + ...)[A:(l -§ + ...)-!] 

When k < 1, 9" ~ O(^), however, for k = 1, 9" ~ 0(^^). For small amplitude 9 (9 ^ 0), 
6'2=^(~ 0(^^^)) is much smaller than 6'fc<i(~ 0(0)). Thus 9' changes much more slowly in the 
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neighbourhood of = for A; = 1. Consequently, the period of oscillation T, for the same 
amplitude is significantly larger. 



TABLE I. Time period of oscillations for different amplitudes at various k values. 



T about = 


k = 


k = 0.5 


k = 0.75 


k = 1 


amplitude vr/lO 


6.32 


8.78 


12.00 


33.71 


amplitude = 7r/20 


6.29 


8.86 


12.41 


66.93 


amplitude = 7r/40 


6.29 


8.88 


12.53 


133.62 


amplitude = 7r/80 


6.28 


8.88 


12.55 


267.12 



Table |T] lists the time periods computed for different amplitudes at various values of k. 
The values in the table reveal a very interesting pattern. For small amplitudes, the time 
period is essentially independent of amplitude when < A; < 1; but for /c = 1, the time 
period increases inversely as the decrease in amplitude. This is explained as follows. For 
< A; < 1, the effective potential is approximately parabolic in the neighbourhood of 6' = 0, 
analogous to the simple pendulum. Thus, the small amplitude motion of the bead is simple 
harmonic, for which the time period is independent of amplitude. However, this parabolic 
approximation does not hold when A; = 1, as the second derivative of the effective potential. 



U (6), vanishes at 6' = 0. Using (12), the time period of oscillation for A; = 1 can be 
calculated as follows, 



T 



de 



^J (cos Q — COS ^^o) (2 — cos Q — cos ^o) 
Expanding the integrand for small gives, 

_ . / dQ 
T 



(17) 



' Ki ~r • • • /V 2! ~r . . . 



2! 4! ' 6! 

Retaining terms up to eighth order in 9q, 
T 



4! 



dO 







+ 



320 (^0 



^0 ^0 



do 



de 



f(i-(|;)^) + i(i-(£n 
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Defining a new variable y 



'o, this can be expressed as, 
dy 



"70 Jo 
8 

^0 Jo 



l-y4 



^0 1 - 

1 - — ^ — — + 



^o' 1 - 



6 1 - ?/4 80 1 - 



-1/2 



Expanding the quantity within brackets in a binomial series and retaining terms upto fourth 
order in 6*0, 



T 



where the constants are given by, 



A 
To 



+ B9o + eel + 0{d, 



10.49 



B 



C 




1 ii-y' 



.&\2 




(19) 
(20) 
(21) 



12(1-7/4)5/2 20 (1 -1/4)3/2 

This shows that for small 6'o, T ^ A/6q, which explains why T doubles when the amplitude 



is halved. Equation (18) gives a very good approximation even if is not close to 0, e.g., 
for 6q = tt/2, the error is only —0.58% (Figure [s]). This can be accounted for by noting that 
the coefficients of higher order terms in 6o decrease rapidly. Thus the expansion converges 
very fast. As an example, for Oq = 1, the third term in (18) contributes only about 0.67% 
of the first term. Another point to be noted is that the time period of oscillation for A; = 0, 
increases with amplitude, whereas, for higher values of k {k > 0.3), the time period decreases 
initially with amplitude (Fig|9]). As Eventually for all k, the time period approaches oo as 
the amplitude approaches vr. 



C. k > 1 



When > 1, in addition to 6 = and 6 = it, a. new extremum of the potential appears 



at fli = cos ^{1/k). From (10), it may be shown. 



[/"(O) = -(A;- 1) < , U"in) = -(k + l) <0 , U" (Vli) 



<:--l>0. 



(22) 
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FIG. 9. Variation of time period with amplitude at various values of k 

6 = n remains a maximum of the effective potential, however, = which was a minimum 
for k < 1, is now a maximum. The minimum is at ^ = fii, where. 



The effective potential is now a double well as shown in Figure [TO] Due to left-right symme- 
try of the hoop, —Qi is also a minimum of the potential. For k > 1, the potential value at 
Qi is less than —1, and the value decreases monotonically as k increases. So, the minimum 
a.t 6 = Qi becomes deeper with increase in k. As k increases from 1 towards oo, fii increases 
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from to 7r/2 in a monotonic manner. The maximum at 6' = is less than the maximum 
at = TT. 

As in the previous section, the motion of the bead in different effective energy regions can 



s > 1 
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FIG. 10. Effective potential U{9) versus 9 iov k> I. 



be identified as follows. 



(i) £ = -1/2(A; + 1/A;) : 

The only two possible positions of the bead are at ±i7i, where both 9' and 9" are zero. 
Here, even though the minima of the effective potential are symmetric; depending on the 
initial conditions, one or the other of these minima will be chosen, leading to a breaking of 
symmetrj3. As the hoop rotates about its vertical axis, the bead simply moves on a circle 
of radius asinfii. 

(ii) -l/2{k + l/k) <e<-l: 

The bead oscillates between 9 = ±a to 9 = ±/3 corresponding to the turning points B, C 
about fii, or 5', C about — This is similar to the oscillations encountered for the cases 
< < 1 and A; = 1, except that, here, the oscillations are confined to one side of the hoop 
(with respect to the vertical axis) and are not symmetric about the minima ±fii. The bead 
swings for longer time below the equilibrium point than above it. As 0' = at a and /3, we 
may write using ([T]), 

cos a + cos (3 = 2 cos Qi or, 

a + (3 a — f3 

cos COS = COS i li 



Hence, cos ^-^^ > cosfii or, i7i — a > /3 — ^li. Thus the turning point B in Figure 
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farther from the point of equihbrium fii, than the other turning point C. The shape of the 



effective potential also confirms this (see Figures 10 and 11). About = 0, the potential is 



symmetric, whereas, about 6 = Qi, it has a steeper ascent towards vr than toward 0. This 



asymmetry in the shape of the potential decreases as k assumes larger values (Figure 11) 



This is illustrated in Figure ^T] where this asymmetry is shown for > 0. In Figure 12 



the variation of the potential energy of the bead as it oscillates about Qi is shown with 
a ~ i7i — 60°. R is the ratio of the time spent below the equilibrium position Qi, to the 
time spent above it. The period of oscillation can be computed as. 



U(e) 




(a) fc = 2, i? = 4.83 (b)A: = 5, i?= 1.58 (c) A: = 20, i? = 1.15 

FIG. 11. Asymmetry in the effective potential about its minimum fii at different k values for 
9>0. 






(a) fc = 2 (b) fc = 5 (c) fc = 20 

FIG. 12. Change in the potential energy of the bead with time at different k values. 



/3 



de 



^ , , ^=^= (23) 

ykJa a/(cos a — cos 6*) (cos ^ — cos /3) 

The trajectory of the bead is confined within a band from 6' = a to 6* = /3, on the surface 



of a sphere of radius a (Figure 13) 



(iii) £ = -1 : 

The motion of the bead is confined within 6' = and a certain f22, where 9' is zero and U 
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(a) T = 25 (b) T = 100 

FIG. 13. Trajectory for k=4, 61(0) = 7r/3, O' (0) = 



matches e (see Figure 10). However, the motion is not oscillatory. If the bead were to start 



at f22, O' {^12) = 0, however, 6" {^12) < 0, as il2 > ^i- As a result, 9' becomes negative and 
hence 9 decreases towards ili. At Qi, 9" = 0, whereas 9' has maximum magnitude but is 
still negative. 



\9: 



\9 



e=ni I 



(24) 



Therefore 9 continues to decrease towards zero. At 9 = 0, 9 = and ^ is zero as well. 
Depending on its initial position, the bead will eventually reach ^ = , either directly, or 
after turning at ±^2, with ever decreasing speeds. Once at this position, the bead will 
remain there until perturbed. The period of oscillation of the bead may be computed by 
direct integration as was done for the case of < A; < 1 with e = 1. However, by invoking 
the time reversal symmetry of the system, one can show that the bead will take an infinite 
amount of time to reach 9 = 0. Let us start with the premise that the bead reaches 9 = 
after a finite time ti. At ti, 9' = 0, and we may equally say that at r = ri, the direction of 
9' has been reversed. Due to time reversal symmetry, the bead must then retrace its path. 
Having reached 9 = however, the bead must stay there forever, as 9' and 9" are both zero 
at this position. Hence the bead cannot reach 6' = in finite time. 



(iv) -1 < £ < 1 : 

The bead oscillates between 9 = ±A corresponding to the turning points D and D' in Figure 



10, where 9 is zero. Here, the bead has enough energy to overcome the local maxima at 



and cross over to the other side. In a single period of oscillation, \9'\ attains its maximum 
value four times, as the bead crosses both the minima fli and —Qi in each half oscillation. 
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The time period can be calculated as 

[' 

/o a/(cos6' - cos A) (2 - A;(cos6' + cos A)) 



4 / . ^ (25) 



(v) 5 = 1 : 

Both^ and^ are zero at 6 = ±7r. Hence, if the bead were at an initial position — tt < 6 < tt, 
it will eventually reach ±tt at ever decreasing speeds. The speed, \6'\ will be maximum at 
±Qi and will have a local minimum at 6' = 0. 



(vi) e>l : 

6' is never zero as e is greater than U{6) for all 6. Thus 6 increases or decreases indefinitely, 
depending on the sign of 9' . The bead will whirl around the entire hoop, in the same 
direction, past its initial position, periodically. The speed of the bead will increase and 
decrease periodically, with a local minimum at ^ = 0, global minimum at ±7r and maximum 
at The trajectory will be similar to the ones shown in Figure [7} 

IV. ENERGY AND CONSTRAINT FORCES 



From ([T]) the total energy of the bead is. 



E = -ma 
2 



a 



(26) 



Recalling the expression for the first integral of motion in ([t]), we may write, 

E = mga ■ e + ma^uj'^ sin^ 9 . (27) 

The rate of change of energy is given by, 

dE 

— = ma^u'^e sin 26 , (28) 
dt ' ^ ^ 

E is an even function of 6. This implies that the net change in energy over a path that takes 
the bead from 6 to —6 or back to 6, is zero, although 6 may not be zero along the path. 
This accounts for the different oscillatory motions of the bead about described in section 

imi 

The constraint forces exerted by the hoop may be calculated using the method of La- 
grange's undetermined multipliers. Introducing two Lagrange multipliers, A and /i, the 
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Lagrangian of the system may be written as, 

L = — [r^ + r^(6'^ + sin^ 6*0^)] + mgr cos 9 + \{r - a) + - cot) , (29) 
with the constraint equations 

r = a and (p = cot . (30) 
Using the Euler-Lagrange equation, the equations of motion for r, and 6 are, 

f = r + 0^ sin^ 6'] + — + cos 6^ , 

m 

d 



mr (psm = /i , 



^2 J, „;„2 

I Y/z/r 

— {r^9) = r'^cj)^ sin 6 cos 9 ~ gr sin ^ . 

Substitution of the constraint equations in the above yield, 

A = —ma [a + co'^ sin^ 9~\ — mg cos 9 , 
/i = ma^uj9 sm.29 . 

A and /i are the constraint force and torque along the radial and azimuthal directions 
respectively^. As the hoop is assumed to be frictionless, it can only exert forces in a plane 
perpendicular to it. Hence, the 9 component of the constraint force must be zero. This 
can be verified by the Newtonian approach. Let Fg denote the constraint force along the 9 
direction. Balancing the forces along 9 direction we get 

Fq — mg sin 9 = ma (9 — sin 9 cos 9) , or, 
Fe = ma ^9 + - sin 9 — co^ sin 9 cos ( 



a 

The factor within the brackets is zero as per (|3|, hence F^ = 0. 
Therefore, the rate of work done by the constraint forces is, 

P = F^{acosm9) = jjiu = ma^co'^9 sm29 . (31) 

which is equal to the rate of change of energy of the bead. 

V. PHASE PORTRAITS AND BIFURCATIONS 

We now study the nature of fixed points and bifurcations, and portray the phase space 
trajectories of the bead- hoop system. By defining a new variable 9i = 9' , may be 
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transformed to a set of two first order equations as, 

9' = 9i, and 9[ = -sm9{l - kcosO) . (32) 

The system described by Q is conservative as there exists a first integral of motion which 
was computed in section |TlJ This conserved quantity represents an energy surface in the 
phase plane. The system motion takes place on this surface at constant height. That is, the 
effective energy remains a constant on the trajectories of the system. Therefore, any isolated 
minimum or maximum of this energy surface is a center and there cannot be any attracting 
fixed points or limit cycled. The system is also reversible, as under the transformations, 

(z) -t, 6^-6, and (ii) t -t, 6^ -6^ (33) 



the governing equations ( p32| remain invariant. This means that if {6{t),6i{t)) is a solution, 
so are {—9{—t), 6i{—t)) and {6{—t), —6i{—t)). The reversibility of the system does not admit 
any repelling fixed points. This can be argued as follows. Let us assume that a repelling 
fixed point does exist. If this is not at the origin, then by virtue of reversibility of the 
system, there must exist an attracting fixed point, located at the reflection of the repelling 
fixed point about 6 or 6i axis. This however, is not possible as the system is conservative. 

Were the origin to be a repeller, then trajectories infinitesimally close to the origin will 
move away from it in all directions. Reversibility requires that for every outward bound 
trajectory close to the origin in the first quadrant, there exists a trajectory moving towards 
the origin in the fourth quadrant. This is in contradiction to the assumption that the origin 
is a repeller. Hence, the fixed points of the bead-hoop system are either centers or saddles. 



Equations (32) show that for all values of k, the fixed points lie on the 9 axis, and in the 
first quadrant, all trajectories move right, as ^' > 0. 

When < A; < 1, there are two fixed points for 6 G [0,7r]. When > 1, an additional 
fixed point appears at ^ = cos~^(l/A;) We first investigate the fixed points and their nature 
by linear stability analysis. The effect of nonlinear terms is considered subsequently. A 



general form of (32) may be written as, 

e' = fie,e,), e[ = gie,e,] 
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A Taylor expansion about a fixed point {6*,6l) yields, 



df 



dO 
dg 



df 



'de 



d9i 
dg 



e = /(r , ei) + + ^-f + 0{e, ri\ er?) + . . . (34) 



'd9i 



V = g{0% ei) + „ + + 0{e, v\ + (35) 



where, ^ = 9 9* and t] = 9i — 9^. For small perturbations about the fixed point {9*,9l), 
quadratic and higher terms in ^ and r] may be neglected. The equations governing ^ and rj 
are given by. 



dl df_ 

ae ae. 



(36) 



/ \ 96* 801 

The matrix on the rhs is the Jacobian matrix at the fixed point {9*, 91). 
The Jacobian matrix at the fixed point (0, 0) is. 



J(0) = . (37) 




As A; — 1 < 0, the eigenvalues A = — k are purely imaginary. Thus (0, 0) is a hnear 

center with angular frequency of revolution given by a; ~ Vl — k. On expanding e in the 
neighbourhood of the origin, 

2e^9l + {l- k)9^ - 1 + 0{9^) , (38) 

it is clear that (0,0) is a minimum of the effective energy. We therefore conclude that (0, 0) 
is a non-linear center as well. 

Thus, the trajectories close to the origin are elliptical with eccentricity e = y/k. The 

trajectories start as circles when k = 0, and become elongated in the 9 direction as k 
increases. For (7r,0), the Jacobian matrix is. 



(39) 



The eigenvalues of this matrix are real with opposite signs, A = ±\/k + 1. Hence, (tt, 0) is 
a saddle with eigenvectors given by, v = (1 A)^. Therefore, (1 -\/F+T)^ is the unstable 
manifold with A = \/k + 1 and (1 — ^/k + 1)^ is the stable manifold with A = —^/k + 1. 
They have equal and opposite slopes, whose value increases from 1 to \/2 as k increases from 
to 1. 
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The phase trajectories for k = 0.75 and k = 1.0 are shown in Figure [14} One can iden- 



tify heterochnic orbit or saddle connection, which is a trajectory connecting two saddles. 
A trajectory starting infinitesimally close to 6 = —it with 6i = 0, crosses the 6i axis at 
a/2(£: + 1). By reversibility, this trajectory must reach 6 = rr, thus forming a heterochnic 
orbit. It may be noted here that 6 = n and 6 = —n correspond to the same position on 
the hoop. The 9i — 6 plane is like the curved surface of a cylinder, with the — vr and vr 



edges joined. The heterochnic orbits in Figure [14| are thus homoclinic orbits on the cylinder 
starting and ending on the saddle at (vr,0). The closed orbits around the center at origin 





(a)fc = 3/4 (b)fc = l 

FIG. 14. Phase portraits for different k values 



represent periodic oscillations about = 0, for — 1 < £ < 1, traditionally called librations. 
The saddle connection or heterochnic orbit represents the delicate motion of the bead where 
it starts near the top of the hoop, swings past the bottom of the hoop and slows to a halt 
as it approaches the top of the hoop again. The trajectories above the saddle connection 
correspond to the whirling motion or complete revolution of the bead around the hoop. 



When A; = 1, the Jacobian at (0,0) becomes, 

(40) 




Consequently, both the eigenvalues vanish. This does not give us much information, so we 
expand the effective energy, e, in the neighbourhood of the origin, 

e=| + ^-l + 0(^^^). (41) 
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The point (0, 0) remains a minimum of the energy surface and hence is still a center, though 
very weak. The weakness of the center is brought out by the elongated shape of the orbits 



in the phase plane (Figure 14). Close to ^ = 0, the decay of 6 is very slow. This lethargic 



decay is termed 'critical slowing down' and sometimes signifies the onset of bifurcatiorP. 



When > 1, a new fixed point appears at (fii, 0). As we noted in section III, the minima 
of U{6) and hence the minima of the energy surface occur at ±Qi. Hence from our previous 
reasoning, (f2i, 0) is a center. This can be further verified by forming the Jacobian at {Qi, 0), 



J(^i) =1 , , • (42) 




J(fii) has purely imaginary eigenvalues A = ±i^yk — 1/k. Thus, (f2i,0) is a center, with 
angular frequency of libration oo ~ ^Jk — Xjk. Let A; = 1 + 5. Then (1 + S)^^ = cosfii. 
Expanding both sides for small S and equating the dominant terms, one finds that Qi ~ 
0{VS). Thus, fii increases much more rapidly than 5 does. 
On the other hand, at (0, 0), 

J(0) = f ° M (43) 
\^ A;- 1 Oy 

now has real eigenvalues A = ±\/k — 1. Therefore, the origin has transformed into a saddle 
having (1 \/k — 1)^ as the unstable manifold and (1 — \Jk — \)^ as the stable manifold. 
They become steeper as k increases. The fixed point at (vr, 0) continues to be a saddle as 
before. 

A trajectory starting at — 0^ and = cuts the Q axes exactly at = ^2- By 
reversibility, we conclude that this forms a homoclinic orbit, enclosing the center at (fii, 0). 
The heteroclinic orbits remain as before. We get closed orbits in the region of the phase- 
plane bounded by the homoclinic and the heteroclinic orbits (Figure 15). By expanding 
the effective energy about (i7i,0), we find that the orbits near the center are approximate 
ellipses described by, 

e\^{k- ^) {e - ^]l)2 = 2^ + (A; - ^) + 0((e - ^^f) (44) 

with eccentricity e ~ ^J\k — 1/k — 1|. These features are shown in the phase portraits in 



Figure 15 The physical interpretation of the trajectories is as follows : the fixed point (fii, 0) 
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9, «• 






(a)fc = l.l (b) fc = 1.1 (central part) (c) fc = 4 

FIG. 15. Phase portraits for k > 1 showing homodinic orbits. 

is the stable equilibrium position. The small orbits around it represent periodic oscillations 
or librations. One can identify homoclinic orbits which begin and end at the same fixed 
point, the saddle at (0,0). This corresponds to the motion of the bead where it slows to a 
halt at 6' = 0. The closed orbits surrounding these correspond to oscillations about = 



with amplitude greater than ^2 of Figure [10] in section |III[ The heteroclinic orbits indicate 
the delicate motion where the bead slows to a halt at the top of the hoop. The trajectories 
outside it represent periodic whirling or complete revolutions of the bead around the hoop. 
These inferences are consistent with the motion of the bead as discussed in section III using 
graphical and differential equation analysis. 

A change in the number and nature of fixed points of a dynamical system, as one of its 
parameters is varied is called bifurcation. When k > 1, the fixed point at (0, 0) changes from 
a center to a saddle. Two new centers emerge on both sides of (0, 0) at ± cos~^(l/fc), moving 
away from the origin as k increases further. Therefore, the system undergoes a supercritical 
pitchfork bifurcation at the origin as k is varied through 1^. Also, |J(0)| = for = 1. 
So, it can be called a zero eigenvalue bifurcation. It is also referred to as a symmetry- 



breaking bifurcation. The bifurcation diagram in Figure 16, encompasses both positive and 



negative values of k. For the bead-hoop system that we have studied, negative values of k 
have no physical meaning. However, there are systems with governing equations similar to 



(32), where negative values of k are allowed. For example, a charged particle moving on a 



vertically rotating circular wire, with a uniform magnetic field Bq in the downward vertical 



direction, has the same governing equations as (32), but with the parameter k given by, 

,2 



(45) 



(g/a) mg 

This allows negative values of > —uj\/uj1, where cu/, is the Larmor frequency. For this 
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FIG. 16. The bifurcation diagram for both positive and negative k. Pitchfork bifurcation takes 
place at A and B. B and B' represent the same point. 



system, negative k values become relevant. In Figure 16, the solid lines denote center and 
the dashed lines denote saddle. As k is made more negative than -1, (±7r, 0) changes from 
a saddle to a center and two symmetrical saddle points fork out towards 9 = ±cos^^(l/A;) 
(pitchfork / symmetry-breaking bifurcation). If we imagine the diagram to be on a cylin- 
drical surface with the top and bottom borders joined end to end, then the figure is very 



symmetrical (Figure 16). 



VI. CONNECTIONS WITH OTHER SYSTEMS 



The system we have studied has interesting connections and similarities with a variety of 
other physical systems. As mentioned previously, the same set of equations are obtained for 
a charged particle moving on a vertically rotating circular wire, in the presence of a uniform 
magnetic field directed vertically downward. However, the parameter k can take negative 
values as well. The phase portraits, nature of fixed points and bifurcations can thus be 
analyzed in a similar manner. 

The system exhibits spontaneous symmetry breaking at = 1 (w = Uc), quite similar 
to a Landau second-order phase transitiorP^. The critical angular velocity of the hoop, oOc, 
is analogous to the critical temperature and the equilibrium position of the bead to the 
order parameter in a Landau system. 
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Another example is the Duffing oscillator described by the equation, 

x + x + ex^ = . (46) 
Expanding (|4]), for < A; < 1, in the neighbourhood of = we get, 

9" + (1 - k)9 +(^-^y^ = 0(9') . (47) 



If we define u = a/I — k r, then the above equation may be written as, 

0+»+6{r^»' = O(''^)- fo^O<..<l. (48) 
Neglecting terms of order 5 and higher in 9, the equation reduces to the Duffing oscillator 



equation (46). Thus the nature of fixed points and bifurcations studied for the bead-hoop 



system in the small amplitude limit, with < k < 1, hold for the Duffing oscillator as 



well. The Duffing equation (46) describes the undamped motion of a unit mass attached 
to a nonlinear spring with restoring force F{x) = —x — ex^. The Duffing equation is also a 
conservative system. The coefficient e, of the nonlinear term is related to the parameter k 
of the bead-hoop system as e = {4k — 1)/(6(1 — k)). As k is varied from to 1, e increases 
from —1/6 to oo. For this range of e, the Duffing oscillator does indeed have a nonlinear 
center at the origirP. 

A variation of our system is the hoop rotating about a horizontal axis. This system can 
operate as a one-dimensional ponderomotive particle trapP. 

The rigid pendulum can be considered a special case of our system with k = 0, and 
many problems in various branches of physics, such as, the theory of solitons, the problem 
of superradiation in quantum optics, and Josephson effects in weak superconductivity can 
be reduced to the differential equation describing the motion of a pendulunP. 

Let us now discuss briefly the case when the hoop is not maintained at constant rotation, 
i,e., is not a constant. Initially, the hoop is given some non-zero and the bead-hoop 
system is left to itself. The total energy of the system is conserved, but both 9 and (p vary 
with time. From the Lagrangian of the system, one can arrive at the following. 



M + m sin^ 9 ' 

ai cos 9 
(1 + 6sin^^ 



9 =-sm9 



(49) 
(50) 
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where M and m denote the masses of the hoop and the bead respectively, and ai = 
kla/{M^g), b = m/M. When 6 < 1, i.e., m<t:M, (p^!^ = u, a. constant. In this 
hmit, the hoop has very large inertia compared to the bead, and therefore a high tendency 
to resist any change in its energy. 



The diverse modes of motion of a bead moving without friction, on a vertically rotating 
circular hoop have been explored using a simple theoretic approach based on symmetry 
arguments and elementary calculus. This simple system exhibits several features of nonlin- 
ear dynamics and hence this study can serve as a good background for investigating more 
complicated nonlinear systems. 
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